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Abstract 

The process of neutrino propagation through an active medium consisting of mag- 
netic field and plasma is analysed. We consider in detail the contribution of a magnetic 
field B into the neutrino self-energy operator S(p). The results for this contribution 
were contradictory in the previous literature. For the conditions of the early universe 
where the background medium consists of a charge-symmetric plasma, the pure B- 
field contribution to the neutrino dispersion relation is proportional to (eB) 2 and thus 
comparable to the contribution of the magnetized plasma. The neutrino self-energy op- 
erator S(p) is calculated also for the case of high-energy neutrinos, which corresponds 
to the crossed field approximation. The probability of the neutrino decay v — > e~W + 
is calculated from the imaginary part of the operator. A simple analytical result is 
obtained for the most interesting region of parameters which was not considered earlier. 
The external magnetic field contribution into the neutrino magnetic moment is calcu- 
lated. The result obtained corrects the formulas existed earlier. We show qualitatively 
that the effect of "neutrino spin light" discussed in the literature has no physical region 
of realization because of the medium influence on photon dispersion. 
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1 Introduction 



The presence of matter or electromagnetic fields modifies the dispersion relation of neutri- 
nos in rather subtle ways because these elusive particles interact only by the weak force. 
However, it was recognized that the feeble matter effect is enough to affect neutrino flavor 
oscillations in dramatic ways because the neutrino mass differences are very small [1,2], 
with practical applications in physics and astrophysics whenever neutrino oscillations are 
important [3]. 

The presence of external fields will lead to additional modifications of the neutrino dis- 
persion relation. There is a natural scale for the field strength that is required to have a 
significant impact on quantum processes, i.e. the critical value 

B e = m 2 /e 1=3 4.41 x 10 13 G . (1.1) 

Note that we use natural units where h = c = 1 and the Lorentz-Heaviside convention where 
a — e 2 /An 1/137 so that e ~ 0.30 > is the elementary charge, taken to be positive. 

There are reasons to expect that fields of such or even larger magnitudes can arise in 
cataclysmic astrophysical events such as supernova explosions or coalescing neutron stars, 
situations where a gigantic neutrino outflow should also be expected. There are two classes 
of stars, i.e. soft gamma-ray repeaters (SGR) [4,5] and anomalous x-ray pulsars (AXP) [6,7] 
that are believed to be remnants of such cataclysms and to be magnetars [8] , neutron stars 
with magnetic fields 10 14 -10 15 G. The possible existence of even larger fields of order 10 16 - 
10 17 G is subject to debate [9-13]. The early universe between the QCD phase transition 
( ~ 10~ 5 s) and the nucleosynthesis epoch ( ~ 10 _2 -10 +2 s) is believed to be yet another 
natural environment where strong magnetic fields and large neutrino densities could exist 
simultaneously [14]. 

The modification of the neutrino dispersion relation in a magnetized astrophysical plasma 
was studied in the previous literature [15-18] . In particular, a charge-symmetric plasma with 
m e <C T <C mw and B <T 2 was considered for the early-universe epoch between the QCD 
phase transition and big-bang nucleosynthesis. Ignoring the neutrino mass, the dispersion 
relation for the electron flavor was found to be [17, 18] 

7tt 2 T 4 / 1 2 \ T 2 eB 



E_ _ V2G F 



15 

(eB) 2 f T 2 



,2 



2ir 2 m 2 w 



In — ^ sin 



,2 



2 — COS0 + 

'u / m w 

(1.2) 



where p is the neutrino momentum and <j) IS the angle between B and p. The first term 
proportional to Gf in Eq. I|1.2J) is the dominating pure plasma contribution [19], whereas 
the second term is caused by the common influence of the plasma and magnetic field [17]. 
The third term is of the second order in (eB/T 2 ) <C 1 but was included because of the large 
logarithmic factor ln(T/m e ) 1 [18]. The dispersion relation of Eq. i|1.2fl applies to both 
u e and v e without sign change in any of the terms. 

The -B-field induced pure vacuum modification of the neutrino dispersion relation was 
assumed to be negligible in these papers. 

However, recently this contribution was calculated for the same conditions [20,21], with 
an absolutely different result: 

^ = v^2G F ^sin 2 ^e-i/( 2 -), (1.3) 
P 8?H 
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where p± is the momentum component perpendicular to the £?-field. It is easy to check 
that this would be the dominant £?-field induced contribution by far and thus would lead to 
important consequences for neutrino physics in media [22,23]. 

Because of importance of the question whether the -B-field contribution into the neu- 
trino dispersion relation was dominating or negligible, an independent calculation of it was 
strongly urged. 

One more promising effect based on using the neutrino dispersion properties in external 
active medium, the so-called "neutrino spin light" , was proposed in the series of papers [24- 
28], however, the medium influence on the photon dispersion was not considered there. 

The paper is organized as follows. We begin in Sees. and [3] with the technique to 
calculate the neutrino self-energy operator by using the charged-lepton, W- and <I>-boson 
propagators in a magnetic field. In Sec. 0] we derive explicit results for the neutrino self- 
energy operator in the limiting cases of a "weak field" eB <C mj and a "moderate field" 
mj <C eB <C Trtyr. In Sec. we find the pure-field correction to the neutrino energy and in 
Sec.[5]we study its possible contribution into the resonance condition for neutrino oscillations 
in the supernova interior. The probability of the neutrino decay v — > e~W + and the neutrino 
magnetic moment in an external electromagnetic field are calculated in Sees. and |H1 In 
Sec. El we study a question whether the effect of "neutrino spin light" has a physical region 
of realization with the photon dispersion in medium taken into account, before concluding 
in Sec. HDJ 

2 Definition of the neutrino self-energy operator 

A literature search reveals that calculations of the neutrino dispersion relation in external 
B- fields have a long history [29-31]. To compare the different results we introduce the 
neutrino self-energy operator S(p) that is defined in terms of the invariant amplitude for 
the neutrino forward scattering on vacuum fluctuations, v — ► by the relation 

M(i/->i/)=-i/(p)£(p)i/(p), (2.1) 

where p is the neutrino four-momentum. Note that we use the signature (+,—,—,—) for 
the four-metric. 

Perturbatively, the matrix element of Eq. I|2.1|l corresponds in the Feynman gauge to 
the Feynman diagrams shown in Fig. 1 where double lines denote exact propagators in the 
external B field. Put another way, the terms in the self-energy operator correspond to these 
Feynman graphs with the external neutrino lines truncated. 

At first glance, the contribution of the diagram with the scalar field should be negligible 
because of the suppression by the factor (me/mw) 2 arising from the coupling of the lepton 
with the scalar $-boson. However, as will be shown later, it is essential in some cases. 

3 Calculation techniques 

The calculation techniques for quantum processes in external electromagnetic fields based 
on exact propagators in the field started from the classical paper by J. Schwinger [32] and 
was developed by A. Nikishov, V. Ritus, A. Shabad, V. Skobelev et al. For a recent review 
see e.g. [33]. 
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Fig. 1. The Feynman diagrams for the v — » v transition in the Feynman gauge. The 
double lines denote exact propagators of the charged lepton, the J-F-boson, and the 
unphysical charged scalar 3>-boson in an external magnetic field. 



The exact propagator for the charged lepton I (for definitcness we take Qi — — e < 0) in 
a constant and uniform magnetic field can be expressed as 



S F (x,y)=j* l {x,y) S (x-y), 
and similarly for the W- and $-bosons: 

G F pa {x,y) = e^y) Gpa{x - y) . 
D F (x,y) = e^'f) D(x - y) . 



(3.1) 

(3.2) 
(3.3) 



where S(x — y), G pa {x — y), D(x — y) are the translationally and gauge invariant parts of 
the propagators. 

The phases Q,(x, y) being identical for all propagators are translationally and gauge non- 
invariant, but they cancel in the two-vertex loop: 



Q{x,y)+Q{y,x) =0. 



(3.4) 



The Fourier transforms of the translationally invariant parts of the propagators are defined 
by: 
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For the Fourier transform S{q) of the translationally invariant part of the lepton propa- 
gator one obtains in the Fock proper-time formalism: 



S(q) 
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where [3 = eB and m,( is the lepton mass, ip a p — F a p/B is the dimensionless field tensor. The 
Lorentz indices of four-vectors and tensors within parentheses are contracted consecutively, 
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e.g. (71/37) = r ) a <-Pai3l^ ■ In the frame where the B field is directed along the 3d axis, four- 
vectors with the indices _L and || belong to the Euclidean {1, 2}-subspace and the Minkowski 
{0, 3}-subspace, correspondingly. For example, p± — (0,pi,p2,0) andpy = (po,0, 0,p 3 ). For 
any four-vectors X and Y we write, both in the invariant form and in the above-mentioned 
frame: 



(XY) l{ = (X<wY)=X Y -X 3 Y 3 , 
(XY) ± = {X W Y) = X 1 Y 1 +X 2 Y 2 , 
(XY) = (XY) l{ -(XY) ± . 

Similarly, for the W- and $-boson propagators in the Feynman gauge we have: 



(3.9) 
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(3.11) 



Magnetic fields existing in Nature probably are always weak compared with the criti- 
cal field for the W-boson, B\y = m^/e ~ 10 24 G. Therefore, the VF-propagator can be 
expanded in powers of (3 as a small parameter: 



G pa {q) 



9pa 



q 2 - rrvjy 



- P 



2<Ppo 



Spa 



(q 2 - m 2 w ) 2 
1 , 2<zi 



(q 2 



(q 2 



' 4 ((pip) 



pa 



(q 2 



m 2 ) 3 



(3.12) 
0((3 3 ). 



Likewise, the asymptotic expression for the lepton propagator S(q) is realised when the 
field strength is the smallest dimensional parameter, j3 <§; m 2 <C m 2 ^. In this "weak field 
approximation" the charged-lepton propagator can be expanded as [34] : 



S(q) = i 



(97) + mi 
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One can see that the contribution of the region of small virtual momenta q 2 ~ mj <C nvfo 



is enhanced in each succeeding term. If the propagator is used for a "moderate field' 
m 2 <S 
taken 



m 2 <C (3 -C rriy^, the expansion is not applicable and the exact propagator (|3.8|) must be 



4 Calculation of the operator 

From the S matrix element for the transition v — > v corresponding to the Feynman diagrams 
of Fig. 1, one extracts the amplitude M.. The neutrino self-energy operator can be written 
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in the form: 
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where g is the electroweak SU(2) coupling constant of the standard model, L = i(l — 75) 
and R = 5(1 + 75) are the left-handed and right-handed projection operators, and 

J$P<P)= J ^S(q)G f3a (q~p), jW(p) = J J0^S(q)D(q-p). (4.2) 

It is convenient to express the structure of the S(p) operator in an external magnetic 
field in terms of the coefficients Al, Bl, Cl, Ar, etc. 



A L {pi) + B L e 2 (pFFj"} + C L e (pF 1 

A R O7) + Br e 2 {pFFj^j + C R e (pFjj 
m v [JCi +i/C 2 e (7F7)] , 



L 
R 



(4.3) 



where F is the external field tensor, and F is its dual. 

Here, the coefficients Ar, Br, Cr, JCi^ are originated from the Feynman diagram of Fig. 
1 with the scalar <i>-boson, while the coefficients Al, Bl, Cl contain the contributions of 
both diagrams. 

The coefficients Al, Ar and ICi, being ultraviolet divergent, do not have independent 
meanings, because they do not give contributions into the real neutrino energy in external 
field at the one-loop level. They are absorbed by the neutrino wave-function and mass renor- 
malization. The coefficients Br, Cr are suppressed by the factor (m„/miy) 2 . The coefficient 
K-2 is suppressed by the factor (mg/mw) 2 , see, however, Sec. OD Thus, the coefficients Bl, 
Cl are of the most interest. 

We present our results for the Bl and Cl coefficients of the £(p) operator (|4.3|) in Table 
1 where the results of previous authors [29-31] and of the above-mentioned papers [20,21] 
are also shown. 



5 Neutrino energy in a magnetic field 

Solving the equation for the neutrino dispersion in a magnetic field (for m v = 0) 



det 



(P7) - B L e 2 (pFF 1 ) L-C L e(pF 7 ) L 



= 0. 



(5.1) 



where the leading terms with Bl, Cl are only included, one obtains for the neutrino energy 
in the field: 
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Table 1. Coefficients in Eq. iA.'M for the neutrino self-energy operator S(p) in an 

external B-field. 



It can be seen that the Bl coefficient gives the main contribution into the neutrino 
energy, because the value C 2 /Bl ~ Gym^ appears to be of the order of the fine-structure 
constant a ~ 1/137, thus leading us beyond the frame of the one-loop approximation. 

Our results strongly disagree with those by E. Elizalde e.a. [20,21]. We think that the 
disagreement arises because these authors use only one lowest Landau level in the charged- 
lepton propagator in the case of moderate field strengths which they call "strong fields." 
However, the contributions of the next Landau levels appear to be of the same order as 
the ground-level contribution [35] because in the integration over the virtual lepton four- 



momentum in the loop the region q 



1 w 



^> (i appears to be essential. 



We confirm the assumption of Refs. [15,17], that the pure magnetic field contribution 
into the neutrino energy does not exceed the plasma contribution. 

For relatively weak field eB <t^. m 2 we find the pure-field correction to the electron 
neutrino energy in a magnetic field and plasma, rewriting Eq. <|1.2H as follows: 
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It is seen from the last terms that the pure magnetic field contribution to the neutrino 
dispersion is proportional to (eB) 2 and thus comparable to the contribution of the magne- 
tized plasma. It is interesting to note that the contributions of plasma and of pure magnetic 
field into Eq. (|5.3|) . containing the electron mass singularities ~ lnm e , exactly cancel each 
other. It looks suspicious that the singularity remains in the case of a moderate field. Pos- 
sibly it means that the plasma term ~ (eB) 2 obtained in Ref. [18] is valid in the weak field 
case only. 



6 Field-induced resonance transition u T}fi — ► v e 

There exists a long-standing problem in the supernova explosion modelling, of searches of 
an energy transfer mechanism from the neutrino outflow to the stalled shock wave for its 
revival [36,37]. 

If a strong magnetic field is generated inside the exploding supernova [9-13], conditions 
could arise for the resonance enhancement of the neutrino oscillations v^ %T — > v e (the effect 
similar to MSW [1,2]), with further i/ e -energy transfer to the stellar matter via the URCA 
processes. 

With the magnetic field contribution into the difference of the neutrino self-energies 
AE = E vt — E Ve (i = (J,,t), when m 2 . <C eB <C m 2 <C m^, the resonance condition for the 
V[ — > v e oscillation is 

2E 3V2n 2 m 2 w \ eB J m N V ' 

where 9 is the mixing angle in the vg, v e system, p is the matter density, Y e is the electron 
fraction with respect to nucleons, is the nucleon mass. It is remarkable that the sign of 
the field-induced neutrino self-energy difference is favorable for the resonance appearance. 

Neglecting the neutrino masses, we obtain the following equation for evaluation of the 
magnetic field strength providing the resonance transition v r — > v e : 

B 2 7 (1 - 0.10 x bxB 17 ) ~ 2.5 x 10 2 x PjJ^l ; (6.2) 

E\o 

where B„ = B/(10 17 G), Pr = p/(10 7 g/cm 3 ), Y . 5 = F e /0.5, E 10 - S/(lQMeV). 

The analysis shows that for realisation of the resonance transition v T>tl — > v e the field 
strength is necessary B > 10 18 G, far exceeding the maximal magnetic field strength which 
is believed to arise inside the exploding supernova. 



7 Neutrino decay v —> e W + in external 
electromagnetic field 

One more interesting result which can be extracted from the neutrino self-energy operator 
is the probability of the neutrino decay v — > e~W + in an external electromagnetic field [30]. 
It is defined by the imaginary part of the amplitude: 

w{v -► e~W + ) = —TmM{v v) . (7.1) 
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In the case of high neutrino energy, which is only interesting for this process, the probability 
is expressed in terms of the dynamical parameter x 



,* = S™, (7.2) 
m w 

and the crossed-field approximation is available. The calculation technique for quantum 
processes in an external crossed field was developed by A. Nikishov and V. Ritus. 

In the paper [30] the probability was written in a general form and in the two limiting 
cases: x 2 ^ A = m e/ m w — 4 1 10 -11 , and x 3> 1. However, rewriting the x parameter in 
the form 

2 _ . o , n -3 f B \ ( E 



X - 3 " 10 UJ [w^ev) ' (7 ' 3) 

one can see that for very wide regions of the magnetic field value and the neutrino energy 
the x parameter belongs to the interval 

A « x 2 « 1 , (7.4) 
which should be much more interesting. 



Our result for the decay probability in the case A <§; x <C 1 is: 



2 G F x 2 
ml 1 + (vs) 



1 - (vs) 



2 mw tan c 



(ts) 



2m w 



(7.5) 



where v = p/E is the neutrino velocity, s is the unit vector of the neutrino spin direction, 
and t = (n x (B x n))/(Bsin</>) is the unit vector lying in the plane of B and n = p/|p|. 
The term with t exists in the case if the neutrino has a transversal polarization. 

The last term in l|7.5|l is provided by the contribution of the Feynman diagram in Fig. 
1 with the charged scalar <E>-boson. It should be noted that the probability is not positively 
defined without it, when the angle between s and n is small, but not equal to zero exactly. 



8 Field contribution into the neutrino magnetic mo- 
ment 



One more case when the contribution of the diagram with the charged scalar <E>-boson is 
essential, is the field contribution into the magnetic moment \i Vl of the neutrino vg. The 
neutrino magnetic moment is expressed in the coefficients of the self-energy operator l|4.3|) 
as follows: 



e m v 



(C L -C R + AIC2) 



■ I) 



In the limiting case x 2 ^ A = mf /m w we obtain: 
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where /iff is the neutrino magnetic moment in vacuum [38,39] 

(0) _ 3eG F m„ 

The leading term of the field correction ~ x 2 m <|8.2[) with a big logarythm coinsides with 
the result of Ref. [30] where the post log terms were not taken into account. The last term 
(1/3) in Eq. I|8.2I) provided by the <&-boson, is relatively small, but it is not parametrically 
suppressed, as it was treated earlier. 



9 No "neutrino spin light" because of photon dispersion 
in medium 

The influence of an active medium on the neutrino dispersion was exploited in the recent 
series of papers by A. Studenikin et al. [24-28], where the so-called effect of "neutrino spin 
light" was discovered. 

The idea was based on the additional Wolfenstein energy [1] acquired by a left-handed 
neutrino in medium: 

E VL ~ E + ^fi- (1 + 4 sin 2 0w) , (9.1) 

E UR ~ Eq, (9.2) 

where N is the number density of background electrons. 

Given the effective vl^rI vertex caused by the neutrino magnetic moment, the decay 
became possible, in those authors' opinion: 

v l * v r ~t~ 7 i (9-3) 

with the photon emission which was called "the neutrino spin light" . 

As is seen from the papers [24-28] , the authors made the kincmatical analysis and calcu- 
lations of the decay probability and other observables, with taking account of the neutrino 
dispersion in medium, but considering the photon created as it was sterile with respect to 
the medium influence, and had the vacuum dispersion, u> = |k|. However, it is well-known, 
that medium modifies essentially the photon dispersion to the form co = |k|/n, where n ^ 1 
is the refractive index. Having in mind possible astrophysical applications, it is worthwhile 
to consider the astrophysical plasma as a medium, where the photon acquires properties of 
a plasmon [40-42]. The dispersion curves for the transversal and longitudinal plasmon are 
depicted in Fig. 2. 

The deviation of the plasmon dispersion in dense matter from the vacuum one is defined 
by the so-called plasmon frequency which is for the relativistic case: 

It should be compared with the Wolfenstein energy defining the neutrino dispersion in 
medium 

G F N /1 , , . 2 \ _ . i c% -it f N 



AEw = i/r (1 + 4 sin2 9w) " L2 eV ( i" 37 -- 3 ' ■ (9 - 5) 
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Here, the scale of the electron number density is taken, which is typical for the interior of a 
neutron star. For smaller densities, the value u p \ exceeds AEw even much greater. 

As can be seen from the dispersion plot, the 4- momentum of the transversal plasmon is 
always timelike, uo 2 > k 2 (n < 1). It means that this plasmon has an effective "mass" which 
is much greater than the energy benefit caused by the neutrino dispersion. So, the decay 
vl — *• VR-ft is kinematically forbidden. This is also true for the decay vl — * ^R"fi in the 
region where the 4-momentum of the longitudinal plasmon is timelike. 

On the other hand, in the region where the 4-momentum of the longitudinal plasmon is 
spacelike, uj 2 < k 2 (n > 1), the decay vl — * vrH is kinematically allowed due to the plasmon 
dispersion (the neutrino Cerenkov process), and the contribution of the neutrino dispersion 
into this effect is negligibly small. It should be mentioned also, that the longitudinal plasmon 
being created in that region, is unstable because of the Landau damping, and the neutrino 
energy is not transformed into the "light" radiation, but in fact into the energy of excitation 
of plasma electrons. Thus, the effect of "neutrino spin light" has no physical region of 
realization. 



10 Conclusions 

• We have calculated the neutrino self-energy operator £(p) in the presence of a magnetic 
field B. In particular, we have considered the weak-field limit eB <C m 2 , and a 
"moderate field" case, m 2 <C eB <C m^. Our results strongly disagree with those by 
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E. Elizalde e.a. [20,21]. We confirm the assumption by J. C. D'Olivo e.a. [15] and by 
P. Elmfors e.a. [17], that the pure magnetic field contribution into the neutrino energy 
does not exceed the plasma contribution. 

• Applying the possible field-induced resonance enhancement of the neutrino oscillations 
to the problem of the supernova shock wave revival, we show that the field strength is 
necessary for this, far exceeding the maximal magnetic field strength which is believed 
to arise inside the exploding supernova. 

• Using the imaginary part of the neutrino self-energy operator, we have calculated the 
probability of the neutrino decay v — > e~W + in an external electromagnetic field. 
We have considered in part the most interesting region of parameters which was not 
analysed earlier: A«^ 2 <1 We have shown that the contribution of the Feynman 
diagram with the charged scalar $-boson is essential, because the probability is not 
positively defined without it. 

• We have calculated the external electromagnetic field contribution into the magnetic 
moment p, Vl of the neutrino vg. We have shown that the contribution of the <I>-boson 
diagram is relatively small, but it is not parametrically suppressed. 

• We have shown qualitatively that the effect of "neutrino spin light" [24-28] has no 
physical region of realization because of the photon dispersion in medium. 

Acknowledgements 

The authors are grateful to G. G. Raffelt and L. A. Vassilevskaya for collaboration and to 
V. A. Rubakov for useful discussion. A. K. is grateful to the Organizers of the XL PNPI 
Winter School on Nuclear and Particle Physics and of the XII St. Petersburg School on 
Theoretical Physics for the invitation to present this lecture and for warm hospitality. 

The work was supported in part by the Russian Foundation for Basic Research under 
the Grant No. 04-02-16253, and by the Council on Grants by the President of Russian 
Federation for the Support of Young Russian Scientists and Leading Scientific Schools of 
Russian Federation under the Grant No. NSh-6376.2006.2. 

References 

[1] L. Wolfenstcin, Phys. Rev. D 17, 2369 (1978). 

[2] S. P. Mikheyev and A. Yu. Smirnov, Yad. Fiz. 42, 1441 (1985) [Sov. J. Nucl. Phys. 42, 
913 (1985)]. 

[3] R. N. Mohapatra and P. B. Pal, Massive neutrinos in physics and astrophysics, Lecture 
Notes in Physics. Vol. 72 (World Scientific, Singapore, 2004), 3rd ed. 

[4] C. Kouveliotou et al., Astrophys. J. 510, L115 (1999). 

[5] K. Hurley ct al., Nature 397, 41 (1999). 

[6] X.-D. Li and E. P. J. van den Hcuvel, Astrophys. J. 513, L45 (1999). 



12 



[7] S. Mereghetti, L. Chiarlone, G. L. Israel and L. Stella, astro-ph/0205122 

[8] R. C. Duncan and C. Thompson, Astrophys. J. 392, L9 (1992). 

[9] G. S. Bisnovatyi-Kogan, Astron. Zh. 47, 813 (1970) [Sov. Astron. 14, 652 (1971)]. 

[10] G. S. Bisnovatyi-Kogan, Physical problems in the theory of stellar evolution (Izdatel'stvo 
Nauka, Moscow, 1989), in Russian. 

[11] S. A. Balbus and J. F. Hawley, Rev. Mod. Phys. 70, 1 (1998). 

[12] S. Akiyama, J. C. Wheeler, D. L. Meier and I. Lichtenstadt, Astrophys. J. 584, 954 
(2003). 

[13] N. V. Ardeljan, G. S. Bisnovatyi-Kogan and S. G. Moiseenko, Mon. Not. R. Astron. 
Soc. 359, 333 (2005). 

[14] D. Grasso and H. R. Rubinstein, Phys. Rept. 348, 163 (2001). 

[15] J. C. D'Olivo, J. F. Nieves and P. B. Pal, Phys. Rev. D 40, 3679 (1989). 

[16] V. B. Semikoz and J. W. F. Vallc, Nucl. Phys. B 425, 651 (1994); 485, 545 (1997). 

[17] P. Elmfors, D. Grasso and G. Raffelt, Nucl. Phys. B 479, 3 (1996). 

[18] A. Erdas, C. W. Kim and T. H. Lee, Phys. Rev. D 58, 085016 (1998). 

[19] D. Notzold and G. Raffelt, Nucl. Phys. B 307, 924 (1988). 

[20] E. Elizalde, E. J. Ferrer and V. de la Incera, Ann. Phys. (N.Y.) 295, 33 (2002). 

[21] E. Elizalde, E. J. Ferrer and V. de la Incera, Phys. Rev. D 70, 043012 (2004). 

[22] E. J. Ferrer and V. de la Incera, AIP Conf. Proc. 689, 63 (2003). 

[23] E. J. Ferrer and V. de la Incera, Int. J. Mod. Phys. A 19, 5385 (2004). 

[24] A. Lobanov and A. Studenikin, Phys. Lett. B 564, 27 (2003). 

[25] A. Studenikin and A. Ternov, Phys. Lett. B 608, 107 (2005). 

[26] A. Grigoriev, A. Studenikin and A. Ternov, Phys. Lett. B 622, 199 (2005). 

[27] A. Grigoriev, A. Studenikin and A. Ternov, Grav. Cosmol. 11, 132 (2005). 

[28] A. E. Lobanov, Phys. Lett. B 619, 136 (2005). 

[29] G. McKeon, Phys. Rev. D 24, 2744 (1981). 

[30] A. V. Borisov, V. Ch. Zhukovskii, A. V. Kurilin and A. I. Ternov, Yad. Fiz. 41, 743 
(1985) [Sov. J. Nucl. Phys. 41, 473 (1985)]. 

[31] A. Erdas and G. Feldman, Nucl. Phys. B 343, 597 (1990). 

[32] J. Schwinger, Phys. Rev. 82, 664 (1951). 



13 



[33] A. V. Kuznetsov and N. V. Mikheev, Electroweak Processes in External Electromagnetic 
Fields (Springer- Vcrlag, New York, 2003). 

[34] T.-K. Chyi et al., Phys. Rev. D 62, 105014 (2000). 

[35] A. V. Kuznetsov, N. V. Mikheev, G. G. Raffelt and L. A. Vassilevskaya, Phys. Rev. D 
73, 023001 (2006). 

[36] J. R. Wilson, in Numerical Astrophysics, edited by J. M. Centrella et al. (Jones and 
Bartlctt, Boston, 1985), p. 422. 

[37] H. A. Bcthe and J. R. Wilson, Astrophys. J. 295, 14 (1985). 

[38] B. W. Lee and R. E. Shrock, Phys. Rev. D 16, 1444 (1977). 

[39] K. Fujikawa and R. E. Shrock, Phys. Rev. Lett. 45, 963 (1980). 

[40] H.A. Weldon, Phys. Rev. D 26, 1394 (1982). 

[41] E. Braaten, Phys. Rev. Lett. 66, 1655 (1991). 

[42] T. Althcrr and P. Salati, Nucl. Phys. B 421, 662 (1994). 



14 



